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We report recent measurements of rare charmless B decays performed by BaBar. The results are based on 
the final BaBar dataset of 424 fb" 1 collected at the PEP-II B-factory based at the SLAC National Accelerator 
Laboratory. 



1. Introduction 

The study of rare B decays is a key ingredient to 
meet two of the main goals of the P-factories: assess- 
ing the validity of the Cabibbo-Kobayashi-Maskawa 
(CKM) picture of CP-violation [l[ by precisely mea- 
suring the elements of the Unitarity Triangle (UT), 
and searching for hints of New Physics (NP), or oth- 
erwise constraining NP scenarios, in processes which 
are suppressed in the Standard Model (SM). In loop 
processes, in particular, NP at some higher energy 
scale may manifest itself in the low energy effective 
theory as new couplings, such as those introduced by 
new very massive virtual particles in the loop Q. In 
NP searches hadronic uncertainties can play a major 
role, expecially for branching fraction measurements. 
Many theoretical uncertainties cancel in ratios of am- 
plitudes, and most NP probes are therefore of this 
kind. 

In the following sections we report recent measure- 
ments, performed by the BaBar Collaboration, that 
are relevant to NP searches in charmless hadronic B 
decays. In Sec. [2] we report the 2008 update of the 
time-dependent analyses of B° — ► r/'K and B° — > 
loKs decays. In Sec. [3] we describe the search for B 
decays to the three-body final states PK$Ks, where 
P = T), T)', 7T°. The Dalitz Plot analysis of B — > Kstttt 
decays is reported in Sec. 2J These measurements are 
related to the CKM angle (3, which is experimentally 
accessible through the interference between the direct 
decay of the B meson to a CP eigenstate and the de- 
cay after B°B mixing. This interference affects the 
time evolution of the decay. 

In the time-dependent analyses of B decays, one of 
the two B mesons produced in e + e~ — > T(45) — > BB 
events is fully reconstructed according to the final 
state / of interest. The flavor and the decay vertex po- 
sition for the other B in the event (B ta g) are identified 
from its decay products. The proper time difference 
between the two B mesons is: 

e -|At|/r 

f(At) = — {1 + q [-rifSf sin(Am d Ai) 

-C f cos{Am d At)}}, (1) 

where r\f is the CP eigenvalue of the final state 
f,q = +!(—!) if the B tag decays as a B° (B ), 



t = (1.536 ± 0.014) ps [|[ is the mean B lifetime, and 
Am d = (0.502 ± 0.007) ps^ 1 is the B° - P° mixing 
frequency [3j. 

The angle can be determined very accurately and 
with small theoretical uncertainties from the analysis 
of the decay rate asymmetries in b — > ccs decays, and 
is related in a simple way to the parameter 

Sec in 

Eq.ru 

Scc = -Vfsin2/3. (2) 

The current world average is 5 CC = 0.67 ± 0.02 0]. 

A parameter S qq — —r\cp sin 2/3 c s (where /3 e ff de- 
notes an effective value of 0) can also be extracted 
from penguin (loop) dominated b — > qqs decays. 
These decays may receive contributions from NP ef- 
fects that can lead to measurable differences AS = 
Sqq — S cc - In some NP scenario, deviations can be 
~ 0(1) [5[. A background to AS measurements as a 
tool for NP searches is represented by the CKM sup- 
pressed tree contributions, which may introduce shifts 
in the parameter S. The SM effect is predicted in most 
models to be a positive shift in AS [6j . The size of this 
shift is related to the ratio of tree to penguin ampli- 
tudes, which depends on the decay mode. Theoretical 
estimates for this ratio are affected by large uncertain- 
ties for most decay modes, with the exception of the 
K S K S K°, (f>K°, and n'K modes. 

The updated measurement of the direct CP- 
violation parameter for B — > K + ir~ is reported in 
Sec.[5j B — » Kit transitions receive contributions from 
tree, color suppressed tree, gluonic penguin (loop), 
and electroweak penguin contributions. Among these 
amplitudes, the dominant contribution is represented 
by the gluonic penguin. The electroweak penguin and 
the color suppressed tree amplitudes are expected to 
be small in the SM, and the direct CP-violation pa- 
rameter for K + ir~ decays should be equal to that 
for K + tt° decays. The rates and direct CP- violating 
asymmetries in the Ktt system reveal puzzling fea- 
tures, that might be indications of NP. 



2. B° -> r)'K° and B° -> uK s 

The BaBar Collaboration has updated the mea- 
surement of CP- violating asymmetries in B decays to 



2 
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rj'K and ujKs using the final dataset corresponding 
to 465 million BB pairs • 

The rj'K channel provides the most precise mea- 
surement of S in a penguin dominated mode. The 
time-dependent analysis of the B — > ?]'K° decays is 
performed on about 2500 signal events, correspond- 
ing to a branching fraction of « 65 x 10~ 6 . The re- 
sulting decay rate distributions are shown in Fig. [TJ 
Both rj'Ks and tj'Kl modes are analyzed, and the re- 
sults are combined in scans of — ln£, taking into ac- 
count the different intrinsic CP parity. With respect 
to the previous measurement, the analysis relies on 
about 20% more data, improved track reconstruction, 
and reconstruction of an additional rf decay channel 
in tj'Kl analysis. The results for the CP-violation 
parameters S and C in this decay mode are 0]: 

S v , K o = 0.57 ±0.08 ±0.02 (3) 
C„'K° = -0.08 ±0.06 ±0.02, (4) 

where the first error is statistical and the second sys- 
tematic. The uncertainties on S and C have decreased 
by about 20% and 25%, respectively, and are still sta- 
tistically limited. The observed discrepancy between 
S cc and S„'ko is not significant. Since this is one of 
the theoretically cleanest modes, the study of decay 
rate asymmetries in B° — > rf K° decays is a bench- 
mark analysis for a Super B factory 0. 

In the B° — ► uiKs channel the tree to penguin ratio 
is not necessarily small. If the two contributions are 
comparable in size they can lead to a non-zero direct 
CP parameter. For the ujKs channel, B daughters are 
reconstructed in the main decay modes, u> — > n + ir~ir° 
and Ks — * ir + n~ . The u mass and angular variables 
are included in the fit. The signal yield for the loKs 
channel is much lower (N = 163) than that for the 
rj'K channel, and the measurement of S and C is 
therefore less precise. The extracted values for these 
parameters are 

S uKa = 0.55t°;*> ± 0.02 (5) 
C uKg = -0.52t°;^±0.03 (6) 

The central value of C is of potential interest as it dif- 
fers from by more than 2a, but more statistics than 
that collected at the B-factories is needed to further 
investigate direct CP-violation in this decay. 
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Figure 1: Data and model projections for n' Ks (top) and 
rj' Kl (bottom) onto At for (a) B° and (b) B tags. Points 
with error bars represent the data: the solid (dotted) line 
displays the total fit function (total signal). In (c) we show 
the raw asymmetry, (N B o — Ngo)/(N B o +N-^o): the solid 
line represents the fit function 
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3. B° -> X°X°P° 



Decays to CP-eigenstatc final states with three CP 
eigenstate particles, two of which are equal, allow 
the measurement of /3 e ff from time-dependent CP- 
asymmetries. Since the CP-eigenvalue is independent 
of the intermediate resonant structure, there is no 
need for an isospin or Dalitz Plot analysis to separate 
the different resonant contributions 9]. BaBar has 
searched for these decay channels in 426 fb _1 data (IoT |. 
A maximum likelihood fit to kinematical and topolog- 
ical variables is performed but no evidence of signal is 
found (see Table |T|, and yields are too small to allow 
a time-dependent analysis. 



Table I Branching fraction results and upper limits at 90% 
confidence level (CL) for B° -> X°X°P channels 0. 



Mode 


B(xlO- G ) 


90% CL UL (xl0~ 6 ) 


7t° KsKs 


2-7lt? ± - 6 


9 


■qKsKs 


2.7j£J ± 1-2 


10 


v'K s Ks 


2.7±i;§ ± 3.4 


20 



4. B° -> K s n 



+ 



IT 



The S parameter can be extracted from a time- 
dependent Dalitz Plot analysis of B — » Kstttt transi- 
tions. BaBar has performed an analysis of this decay 
using a data sample of 383 million BB pairs fill ]. 

The decay amplitude can be parameterized as 
a function of two Mandelstam-likc variables s+ = 

according to the isobar 



model: 



and 
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A{ 8+ ,eJ) = ^| Cj |e-^ x (7) 

3=1 

Ri (m)X L (\f \r')X L (\q\r)T 3 (L,p, q). 

Here the complex coefhcients \cj \ and e~ l ^ j are the rel- 
ative magnitude and phase for the channel j, Rj(m) 
is the lineshape term, Xl are Blatt-Weisskopf bar- 
rier factors [12j . Tj is the angular distribution, p 
is the momentum of the bachelor particle, q is the 
momentum of one of the resonance daughters, L is 
the orbital angular momentum between p and the 
resonance momentum, and a star denotes the B- 
meson rest frame. The signal model includes both 
tttt (p(770),/o(980),/a(1270),/ x (1300),xoo) and Kn 
(K* (892), (Kn)s- wave) resonances as well as non- 
resonant contributions, for a total of 15 complex scalar 
amplitudes. 

The choice of s+ and s_ as independent variables 
is impractical because both the signal events and the 



Figure 2: Standard (top) and square (bottom) Dalitz plots 
of 5° -> K S tt + -k~ candidates [ll[. 



combinatorial e + e~ — > qq events populate the kine- 
matic boundaries of the Dalitz Plot. The following 
new coordinates are therefore introduced: 



1 

7T 



m = — arccos 2 



- 1 



1 



-co, 

7T 



where mo is the tt + it~ invariant mass, to™ 3 " 35 and to™ 111 
are the kinematic limits of m , and is the 7r + 7r _ res- 
onance helicity angle. The following transformation is 
applied, 



ds+ds- 



I det Jldm'dO' 



(8) 



which defines the square Dalitz Plot shown in Fig. [51 
The distribution of the proper-time difference be- 
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Figure 3: Distributions of At when the Bt ag is a B° (top), B° (middle), and the derived At asymmetry (bottom). Plots 
on the left (right) hand side, correspond to events in the f o (980) Ks (p(770)° Ks) region [ill ]. 



tween the two _B-mesons is given by: 

e -|A*IA _ 

f(At) = [L4| 2 + L4| 2 

T 

=F {\A\ 2 - \A\ 2 ) cos(Am d At) 
±77/21111 [AA*] sm(Am d At)] , 

where the CP- violation parameters are expressed in 
terms of the amplitudes A (for B decays) and A (for 
B decays). From the results of this analysis for B° 
decays to p(770)°K s and f (980)K s CP-eigenstates 
it is possible to calculate /3 e ff as: 



0eS = ^ arg (c k c* k ) 



(9) 



The direct and mixing- induced CP asymmetries can 
be calculated as: 



Cl 



s k 



cfe 



Cfc 



|Cfe| z + |c fc | 
21m (c* k Ck) 



Cfe 



|Cfc 



(10) 
(11) 



Signal enriched distributions for At and At- 
asymmetry for events in the regions of yo(770) and 
/o(980) are shown in Fig. [3] Two solutions for the 
S and C parameters are present, almost degenerate 
in likelihood [ll|. For the /o(980)i^s channel, there 



is an evidence of non-zero <S* at about 3.5cr signifi- 
cance (see Fig. |4]). Including systematic and Dalitz 
plot model uncertainties, the confidence intervals at 
95% confidence level (CL) for the measured value of 
/3 off in the B° -> p°(77Q)K s and B° -> 70(980)^5 
channels are [TT[ | : 

- 9° < (3 cS (p (77Q)K s ) < 57°, (12) 
18° < Mfo(980)K s ) <76° (13) 

These measurements are consistent with the SM. 

The analysis of B° — * Kstt + tt~ decays provides 
further non-trivial constraints on the p — 7/ plane 
(CKM angle 7). In particular, the phase difference 

between isobar amplitudes for B° and B decays to 
K*(892)+n- and (Kir)+_ wave 7r- 



A<t>K*TT 



argc x . (5 
arg c. 



Kn)tir- 



- arg 0^.(892)-^+ 

ar gC( i f 7r )- 7r + 



can be used to extract information about 7 1 1 31 ] - The 
one-dimensional scans for AcpK'n and Atp^^s* are 
presented in Fig. [5J Due to the small overlap between 

the phase space regions accessible to B° and B , the 
sensitivity to these phases is limited. The available 
statistics allows to exclude only small regions at 95% 
CL. Including systematic effects, the excluded regions 
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Figure 4: Two-dimensional scans of — 2A log £ as a function of (5, C), for the / (980)if s (left) and p(770)°K s (right) 
isobar components The value — 2 A log £ is computed including systematic uncertainties. Shaded areas, from the 

darkest to the lightest, represent the one to five standard deviations contours. The • (*) marks the expectation ased 
on the current world average from b — > ccs modes [J] (zero point). The dashed circle represents the physical border 
S 2 + C 2 = 1. 




Figure 5: Statistical (dashed line) and total (solid line) scans of — 2Alog£ as a function of the relative phases 
A$(A'*(892)tt) (left) and A$((Kn)o) (right) [llj. Horizontal dotted lines mark the one and two standard deviation 
levels. 



are Oil 



5. B c 



- 137° < A<^-.( 89 2 )7r < -5°, 
-132° < A0 (/f7r)s7r <+25° 

K+7T- 



The direct CP violation Aqp parameter for B — > 
K + tt~ has been extracted from a fit to the 7r + 7r~, 
K^tt^, and K + K~ final states, using 467 million BB 
pairs 14]. Since K + ir~ is a self-tagging mode, Aqp 



is obtained by simple event-counting: 

K-TT+) - N(lf -> K+TT-) 



A C p 



N(B* 



N0 3 



K-TT^ 



N(Tf 



K+TT-) 



(14) 



This analysis yields a 6.1a observation of direct CP 
violation, A CP {K+^-) = -0.107 ± 0.016l^oo4 H- 
By combining this result with Ac<p(K + ir°) and taking 
the averages with Belle [J], one obtains: 



A CP (K+n-) 
A CP (K + n°) 



= -0.098jffl 
= +0.050 ±0.025. 



(15) 
(16) 



G 
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This result corresponds to a 5er deviation from the SM 
expectation. 

Several interpretations have been proposed for this 
discrepancy. In some models this result is considered 
as a hint of NP, which might manifest itself with large 
contributions from electroweak penguins [l5j]. Stan- 
dard Model explanations have also been formulated, 
which involve lar ge c olor-suppressed trees and nonper- 
turbative effects 11611 . 



6. Conclusions 

In this paper we have reported recent results on 
charmless hadronic B decays obtained by the BaBar 
experiment. These results can be used to set non- 
trivial constraints on the Unitarity Triangle, most no- 
tably on sin 2/3 and 7. 

Penguin dominated decays provide several observ- 
ables that can be used as NP probes. No hint of 
NP has been found in the channels described in this 
report, most deviations being within 2a from Stan- 
dard Model predictions. For the theoretically cleanest 
modes, the experimental uncertainty on AS measure- 
ments is dominant, thus demanding for more statistics 
and precision. 

In B — > Kit decays a 5tr discrepancy with respect to 
naive Standard Model expectations is observed. The 
interpretations range from SM explanations to NP ef- 
fects. 
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